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In this paper, we present results of extensive Monte Carlo simulations of diffusion- limited aggre- 
gation (DLA) with a seed placed on an attractive plane as a simple model in connection with the 
electrical double layers. We compute the fractal dimension of the aggregated patterns as a func- 
tion of the attraction strength a. For the patterns grown in both two and three dimensions, the 
fractal dimension shows a significant dependence on the attraction strength for small values of a, 
and approaches to that of the ordinary two-dimensional (2D) DLA in the limit of large a. For non- 
attracting case with a = 1, our results in three dimensions reproduce the patterns of 3D ordinary 
DLA, while in two dimensions our model leads to formation of a compact cluster with dimension 
two. For intermediate a, the 3D clusters have quasi-2D structure with a fractal dimension very 
close to that of the ordinary 2D-DLA. This allows one to control morphology of a growing cluster 
by tuning a single external parameter a. 

PACS numbers: 61.43.Hv, 47.57.eb, 66.10.C-, 68.43.Jk, 82.40.Ck 



I. INTRODUCTION 

By immersing an object with large surface-area-to- 
volume ratio into an ionic solution, it is surrounded by a 
double layer of electrical charge that significantly influ- 
ences its surface behavior [1, 2]. Electrical attraction of 
free ions with thermal motion in the fluid by the surface 
charge, forms a second layer, known as diffuse layer, that 
shields out the Coulomb potential of the surface layer 
and makes the whole structure electrically neutral. Ap- 
plications of double layer range from plasma physics [3], 
to colloidal science [4], and micro- and nano-fluidics [5]. 

As a simple theoretical model for the formation of 
a diffuse layer in a special case, we implement exten- 
sive Monte Carlo simulations for attracted random walks 
(ARW) introduced by Saberi [6] , in order to study the for- 
mation of ionic aggregates on and near an infinite plane 
located at z = 0. An infinite surface charge exerts a 
uniform constant electric force on the free ionic particles 
in the fluid. It is, therefore, reasonable to consider an 
infinite attractive plane of uniform attraction strength a 
which acts on a free random walker launched from a point 
far from the attractive plane. An ionic seed is located at 
the center of the infinite plane. Upon contacting, the 
free ionic random walker sticks irreversibly to the cluster 
due to an absorbing bond introduced between the ionic 
particles. 

The diffusion of an attracted Brownian particle of mass 
m in a fluid may be given by the following Langevin equa- 



tion 
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where the second term denotes a viscous-like friction 
force with drag 7, the third term stands for an attraction 
force of strength a exerted by an infinite plane at 2: = 0, 
and ^{t) is a Gaussian white noise characterized by 
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with D being the diffusion coefficient. 
While such a Langevin-like equation (1), to our best 
knowledge, has not been studied previously, we imple- 
ment such a process to produce a diffusion-hmited ag- 
gregated pattern on an attractive plane and investigate 
its fractal properties in terms of the strength of attrac- 
tion a. We apply two different mechanisms in our model 
to produce the aggregates. First, the aggregates are let 
to grow freely in three dimensions, i.e., the 3D- ARW 
sticks to the aggregate upon touching the cluster giving 
rise to the formation of a 3D pattern. Second, the ag- 
gregates are restricted to only grow within the attractive 
plane, i.e., the 3D- ARW sticks to the aggregate only if 
it touches the cluster within the attractive plane, which 
leads to the formation of a 2D pattern. 
Our results indicate that for small a, the fractal dimen- 
sions for both cases depend significantly on the strength 
of attraction, and in the limit of large a, they converge 
to the fractal dimension of the ordinary two-dimensional 
(2D) diffusion-limited aggregation (DLA). This rapid 
convergence, however, implies the formation of a quasi- 
2D pattern for a 3D aggregate for intermediate vales of a. 
Furthermore, for the two limiting cases, i.e., a = 1 and 
a ^ 00, our model reproduces the results of the ordinary 
DLA. Loosely speaking, for a = 1, the 3D- ARW is identi- 
cal to the ordinary 3D-RW and thus our model generates 
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the same patterns as the ordinary DLA in three dimen- 
sions. In the hmit of a oo, the 3D-ARW can only 
move within the attractive plane and is indeed identical 
to the 2D-RW and thus, in both cases, the patterns of 
the ordinary DLA in two dimensions are reproduced. The 
only difference is for a = 1 in two dimensions, where our 
model gives a new class of aggregates which are compact 
clusters with dimension two. 

The ordinary DLA, known also as a Brownian tree, was 
originally introduced by Witten and Sander [7] to model 
the aggregates of metal particles formed by adhesive con- 
tact in low concentration limit. This model is shown to 
describe many pattern formation processes including di- 
electric breakdown [8], electrochemical deposition [9, 10], 
viscous fingering and Laplacian growth [11]. According 
to this model, the walker does a random walk process ini- 
tiated from infinity like as described above in our model 
but without experiencing any external force. Many of its 
fractal properties have then been reported from various 
simulations done both in real and mathematical plane 
(see for example [12] and references therein). 

This study uncovers a new insight into various related 
phenomena with a discrete time lattice walk [13, 14], 
relaxation phenomena [15], exciton trapping [16] and 
diffusion-limited reactions [14, 17]. 



II. ATTRACTED DIFFUSION-LIMITED 
AGGREGATION (ADL A) -DETAILS OF 
SIMULATION 

In this section, we investigate aggregation of ARWs 
on and near an attractive plane based on the algorithm 
proposed in [7] for ordinary DLA clusters, accompanied 
by the details of simulations. For M being the total 
number of particles, or equivalent ly the cluster mass^ 
the simulation box is considered to be a cube of size 
= Ly = 20VM and = 60VM. The box is di- 
vided into a regular lattice of Lx x Ly x L^ cubic cells of 
unit size. Periodic boundary conditions are applied along 
all three directions. The attractive plane is a horizontal 
plane at z = ^ := whose bottom- left corner is taken as 
origin of the coordinate system. We store information of 
each lattice site by an integer in the computer. Accord- 
ingly, the total memory to be occupied by the lattice will 
be equal to ALxLyLz ^ lO^M^/^ Bytes. One sees that 
as the total number of particles M increases, the total 
memory grows rapidly and for large M the total memory 
becomes so huge that can not be handled by conventional 
computers. In order to prevent excessive memory usage, 
we apply different restrictions to our algorithm to reduce 
the memory. These tricks will separately be discussed 
for both 2D and 3D clusters in Sec. II A and Sec. II B 
respectively. 

First, let us briefiy review the probabilities governing 
the movement of an ARW [6] . For the random walker at 
z 7^ 0, there are six directions which can be chosen by the 
walker. One is that points towards the plane whose prob- 



ability is set to ap, where a > 1 is strength of the attrac- 
tion. For the five directions being left, the correspond- 
ing probability of each link is equally set to p, providing 
p = For the random walker at z = 0, there are two 
directions perpendicular to the plane whose probability 
is p' ^ and four directions on the plane each of which with 
the probability ap' ^ such that p' = j^^- This setting of 
probabilities guaranties that the walker likes to ramble 
on and around the attractive plane. 

It is well known that formation of DLAs is domi- 
nated by diffusive motion rather than the convective one. 
Therefore, the walker must start to move far enough away 
from the attractive plane. Once a seed particle is set- 
tled on the attractive plane, a random walker is intro- 
duced on a random lateral position at height zq = 15a/M, 
and moves according to the probabilities having been de- 




FIG. 1: Snapshots of the growing clusters for a = 1.2 (Top), 
and a = 10 (Bottom), with the same number of particles 
M = 10^. The radius of gyration is Rg — 251.5 and Rg — 
410.54 for the top and bottom figures, respectively. Distance 
from camera is the same for both clusters. One sees that 
the gyration radius of the clusters, Rg, increases with the 
attraction strength a. Moreover, for small a the cluster has 
fat arms, while for larger a the chains are much elongated and 
delicate. The pictures are rendered by the Povray software. 
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scribed above. The walker keeps moving until it sticks 
to the seed. Following that, a new random walker is in- 
troduced at zq = 15^/M from a random lateral position, 
and is let to move until it sticks to one of those par- 
ticles being frozen. This procedure is applied until the 
desired number of the particles in the cluster, i.e., M, is 
reached. For the growth of 2D and 3D structures, we im- 
plement different additional rules that will be discussed 
in Sec. II A and Sec. II B. 



A. 2D structures 

2D structures are formed based on a land- escape pro- 
cedure. When a particle is either atz = +lorz = — 1, 
it may land on the plane if its shadow on the plane is 
not occupied, otherwise it must escape to one of the five 
remaining directions with equal probabilities. Once the 
particle has landed on the plane, it will diff"use accord- 
ing to the aforementioned rules until it finds at least one 
frozen particle on its neighborhood and gets stuck to that 
point. It should be stressed out that according to the 
probabilities governing the motion of the ARWs, there 
still exists a probability for the walker to break from the 
plane and fiy back into the space. 

From our numerical experiments, we empirically found 
that the radial distance of the farthest particles from the 
seed never exceeds r = S^/M limit. We used this fact to 
accelerate our simulation runs, in a way that the attrac- 
tive plane is being divided into two distinct regions, the 
inner region where the cluster forms, and, the outer re- 
gion where the random walker moves freely with no need 
for exploring a possible occupied nearest-neighbor site. 
We employed this fact in our simulations and found that 
it considerably decreases the computational effort. 

In figure 1, two snapshots of the growing clusters on 
the attractive plane are shown for two different strength 
of attraction a = 1.2 and a = 10, shown at the top 
and bottom of the figure, respectively. The cluster mass 
M = 10^ and the distance from camera are considered 
to be the same for both patterns for a comparison. 
Both figures have a fractal structure with a rotational 
symmetry about the 2;- axis and the same number of main 
branches. For lower a, the particles are being condensed 
along the main arms giving rise to much fatter branches 
than for the larger a for which, the cluster has deli- 
cate chains elongated on the attractive plane with deeper 
fjords and sharper tips. It is also evident that the gyra- 
tion radius Rq of the clusters increases with the attrac- 
tion strength a. Here, the gyration radius —in units of 
a particle diameter— is found to be Rg = 220.8 and 
Rg = 410.54 for a = 1.2 and a = 10, respectively. 

In order to get a more quantitative insight about the 
structure of the clusters, we compute the gyration radius 
Rg of different clusters grown under the same conditions 
but with different strength of attraction a. Figure 2 il- 
lustrates the cluster mass M as a function of the gyration 
radius Rg of the clusters for various a. Since the ensem- 
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FIG. 2: Cluster mass M as a function of the radius of gyration 
Rg of the 2D growing clusters on the attractive plane for 
different values of ex. For a > 1.15, the curves exhibit a full 
power- law behavior Eq. (3). 




FIG. 3: A typical example of a 2D cluster grown on the attrac- 
tive plane with a = 1. The number of particles is M = 1875 
which forms a compact cluster of dimension ?^2, with gyration 
radius Rg ^ 19.62. 



ble averaging plays a crucial role in such computations, 
the averages, for each point in the Fig. 2, are taken over 
90 different clusters for given values of a and M. It is 
well-known that the following relation is usually held for 
fractal patterns 



M (xR^^, 



(3) 



where Df is the fractal dimension of the cluster. Our 
results are in accord with this power-law relation (3) for 
a > 1.15. For smaller values of a, the ARW crosses over 
to the ordinary 3D-RW and starts to fill in the empty 
regions between the branches. This leads the cluster to 
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being much condensed. A typical example of such a com- 
pact cluster for a = 1 is shown in Fig. 3. Although the 
power-law behavior is disturbed, within our numerical ac- 
curacy, for a — )► 1, their outer perimeter may keep their 
fractal properties as it comes out from Fig. 3 [22]. This 
would be investigated in more detail in a future work. 

According to our model, for a = 1 the attraction 
strength is set to zero. As a principle of DLA, random 
walkers should initiate their motion from the infinity (far 
from the seed). Therefore, there is a very little chance to 
find a random walker, triggered at infinity, on the surface 
of the attractive plane. Within this small chance, again 
there is a very little probability to keep the particle mov- 
ing on and around the plane. Once the particle lands on 
the plane, normally it tends to escapes to infinity and 
hardly finds a frozen particle in its neighborhood to join 
to the cluster. Consequently, for 2D clusters at a = 1 the 
growth rate is almost zero and we observe only clusters 
with small mass 500 < M < 3000 which is not enough to 
examine any scaling relation. Therefore, we get a poor 
statistics for a = 1 (see also Fig. 6). 



B. 3D structures 

In contrast to 2D structures, the growth model for 
3D clusters is not restricted to the surface of the attrac- 
tive plane. Instead, the random walker in this case can 
get stuck in any lattice site whose each of the nearest- 
neighbor sites is already occupied, whether or not being 
in the attractive plane. 

In figure 4, we show some of the resulting patterns 
obtained by this algorithm in three dimensions. The 
figure shows a side view of the grown clusters consist- 
ing M = 10^ particles with a = 1.2,2 and 10, with 
Rg = 199.6,342.8 and 405.9, from top to bottom, re- 
spectively. Two different colors are chosen for the par- 
ticles: the metallic (gray) color for those sitting on the 
attractive plane at z = 0, and the red (dark) color for 
particles sitting elsewhere with z 0. As can be seen 
from the figure, the number of red (dark) particles is sig- 
nificantly dependent on the strength of attraction. For 
small a, the ARW can easily wander around and thus, 
it is more likely for the particle to get a chance to stick 
to the cluster somewhere at z 7^ 0, while for larger a, 
the ARW is almost rambling on and near the attractive 
plane which increases the sticking probability within the 
plane. In the limit of a 00, the 3D- ARW falls onto 
the ordinary 2D-RW giving rise to the formation of ordi- 
nary DLA clusters in pure two dimensions. For a = 1, at 
the other limiting side, the 3D-ARW is the same as the 
ordinary 3D-RW and we get back the ordinary 3D-DLA 
patterns. 

We empirically found that for 1 < a < 1.2, the thick- 
ness of 3D clusters on either side of the attractive plane 
never exceeds 35 particle diameter, even for clusters of 
size up to M = 10^. Also, for a > 1.2, the height of 
the tallest tower in the cluster never reaches 10. Based 




FIG. 4: Three typical examples of 3D grown clusters of num- 
ber of particles M = 10^. The metallic (gray) particles lie 
on the attractive plane, and the rest— not on the surface— are 
illustrated by the red (dark) color. The strength of attraction 
is considered to be a = 1.2, 2 and 10, with Rg = 199.6, 342.8 
and 405.9, from top to bottom, respectively. The maximum 
height of the towers perpendicular to the attractive plane 
never exceeds 10 particle diameter on either side of the at- 
tractive plane, implying the formation of quasi-2D clusters. 



on this experimental evidence, similar to that of the 2D 
case, the simulation box can technically be divided into 
two following distinct volumes to reduce the memory. An 
inner volume i.e., with —35 < z < 35 for 1 < a < 1.2, 
and — 10 < 2: < 10 for a > 1.2, where the cluster can 
form and the rest as the outer volume where the particle 
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FIG. 5: Cluster mass M as a function of the radius of gyration 
Rg of the 3D growing clusters formed around the seed placed 
on the attractive plane for different values of a. 



In figure 6, we compile our data for the fractal dimen- 
sions of 2D and 3D clusters which are represented by the 
open circles (Q) and squares (□), respectively. 
For both 2D and 3D clusters, the fractal dimension sig- 
nificantly depends on the strength of attraction a, es- 
pecially for a < 2. The smaller the a, the denser the 
clusters are. For larger a, the fractal dimensions in both 
cases, are less dependent on a, rapidly converging to the 
value Df ^ 1.72 which, within the statistical errors, is 
almost in accord with that of the ordinary DLA clusters 
in two dimensions with Df ^ 1.71. 

As discussed before, the overall intermediate fractal 
behavior seems to be dominated by the two crossover lim- 
its in the statistical behavior of the underlying diffusion 
process i.e., the crossover of the 3D-ARW to ordinary 
3D-RW and 2D-RW for a = 1 and a ^ oo, respectively. 
For a = 1, the 2D clusters have a compact structure of 
dimension ^ 2, and for the 3D clusters we obtain a frac- 
tal dimension Df ^ 2.53 very close to the expected value 
for the ordinary 3D-DLA [19]. 



moves freely without an extra check for looking up for an 
occupied nearest-neighbor site. We applied this fact in 
the structure of our programming codes and found that 
it considerably speeds up the simulations. 

The quantity of interest is again the gyration radius, 
Rg, for clusters of different mass M. Figure 5 shows 
our results of cluster mass for 3D clusters as a function 
of their gyration radius for different attraction strength 
a. Each point in Fig. 5 is averaged over an ensemble of 
90 independently grown samples. The mass of clusters 
range between 500 < M < 32000. Error-bars for both 
Rg and M are visible in the plot. We found a perfect 
scaling behavior for a > 2, whereas for 1 < a < 2 our 
data fits to a power-law with higher uncertainty (this is 
also much evident from the inset of Fig. 6). One also 
sees that for a constant cluster mass, the gyration radius 
increases upon increasing the attraction strength which 
is due to the elongation of the growing cluster within the 
attractive plane. As mentioned above, the thickness d 
of the clusters along the 2;-axis never exceeds 70 particle 
diameter even for clusters of size up to M = 10^. This 
implies a considerably small relative thickness, defined as 
d/RG <C 1, for the patterns, and thus indication of the 
formation of quasi-2'D clusters. This would later be con- 
firmed in the next section where we compute the fractal 
dimension of the grown patterns. 



III. FRACTAL DIMENSION OF 2D AND 3D 
ADLA CLUSTERS 

In this short section, we compute the fractal dimension 
of the clusters as a measure characterizing the statistical 
complexity of the grown fractal patterns in both two and 
three dimensions. The fractal dimension of the clusters 
can easily be calculated from the data given in Fig. 2 and 
Fig. 5 by examining the scaling relation (3). 



IV. DISCUSSION AND CONCLUSION 

In conclusion, we introduced a model of aggregation 
based on an extension of diffusion-limited aggregation 
where the underlying diffusion process is a 3D Brownian 
motion (or, a random walk on the lattice) which is at- 
tracted by a plane with strength a. The seed particle is 
placed on the attractive plane. 

In two dimensions, the fractal properties of the aggre- 
gated cluster is shown to be dependent on the attrac- 
tion strength giving rise to formation of patterns with 
1.71 < Df < 2. The two limiting values for the fractal 
dimension Df/m the scahng region, are discussed to be 
governed by the crossover of the 3D- ARW to the ordinary 
3D-RW for a ^ 1 from one side, and its convergence to 
the ordinary 2D-RW at the other side with a ^ oo, which 
leads to ordinary 2D-DLA patterns of Df 1.71. 

In three dimensions, the rotational symmetry, which 
is present in the ordinary 3D-DLA clusters with a = 1, 
is broken by the attractive plane with a > 1, and the 
model leads to formation of clusters with 1.71 < Df < 
2.53. For intermediate a, we obtain quasi-2'D clusters 
whose relative thickness, defined as the thickness of a 
cluster perpendicular to the attractive plane rescaled by 
its radius of gyration, is relatively small. Our results 
indicate a scaling region with a > 2, in which the fractal 
structure of the clusters can be characterized by a single 
fractal dimension Df ^ 1.72, independent of a. 

Recently, there has been a growing attention and 
frequent reports on the fact that morphology of self- 
assemblies can be controlled by tuning various physical 
parameters such as the packing fraction (/), and the attrac- 
tion strength. This is specifically important for handling 
of soft materials and food processing. In our study, we 
found out that a crossover from ramified, fractal clusters 
to compact aggregates occur upon decreasing the attrac- 
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FIG. 6: The fractal dimension Dfoi the attracted diffusion- 
hmited aggregation clusters grown in two and three dimen- 
sions, represented by open circles (Q) and squares (□), re- 
spectively, as a function of the attraction strength a. The 
inset shows the corresponding error bars within which the 
power-law behavior Eq. (3) holds. 



tion strength a at the limit of small attraction strength. 
For the two dimensional clusters the threshold is found 
to be ac = 1.15. As a result, in our model morphologies 
of the growing clusters can be controlled by tuning the 
attraction strength a. 

In a recent study of aggregation of hard spheres [20] , 
a similar crossover from a compact to ramified aggre- 
gation is found which in the former regime the scaling 
relation fails while in the latter case the scaling holds. 
The crossover from compact to fractal clusters happens 
upon increasing the packing fraction (j). The threshold 
is found to be at (j) — 0.55 which is very close to the 



glass transition point at (\)g = 0.58. Therefore, in this 
case, the structure of aggregations can be controlled by 
the packing fraction 0. Zhang et.aL, [21] report results 
of experiments on colloidal suspensions with short-range 
attraction and long-range repulsion. They show that the 
morphology of the clusters can be controlled by tuning 
either the attraction strength or the packing fraction cj) 
resulting to an elongated- to-branched crossover. 

These reports imply that the compact-to-ramified 
crossover may be considered as a universal property of 
aggregation of particles. Consequently, this can suggest 
the possibility of proposing a single phase-diagram that 
captures such behavior for different systems with differ- 
ent underlying physical processes and inter-particle in- 
teractions. 

A possible future work for constructing a much more 
realistic model of diffusive layer may be an extension of 
the present model by considering a number of mobile 
seed particles on the attractive plane that the growth 
of each can be given by our implemented procedure. In 
that case, the growing clusters can stick together upon 
contacting. The results will appear elsewhere in the near 
future. 
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